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ABSTRACT 



Aims. We explore the gamma-ray bursts' (GRBs') contributions in constraining the dark energy equation of state (EOS) at high 

(1.8 < z < 7) and at middle redshifts (0.5 < z < 1.8) and estimate how many GRBs are needed to get substantial constraints at high 

redshifts. 

Methods. We estimate the constraints with mock GRBs and mock type la supernovae (SNe la) for comparisons. When constraining 

the dark energy EOS in a certain redshift range, we allow the dark energy EOS parameter to vary only in that redshift bin and fix EOS 

parameters elsewhere to -1. 

Results. We find that it is difficult to constrain the dark energy EOS beyond the redshifts of SNe la with GRBs unless some new 

luminosity relations for GRBs with smaller scatters are discovered. However, at middle redshifts, GRBs have comparable contributions 

with SNe la in constraining the dark energy EOS. 
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1. Introduction 



Since the discovery of t he accelerating expans ion of the uni- 
verse ('R iess et alJll998l: IPerimutter et al1ll999 '). a lot of work 
has been done on constraining the behavior of dark energy. In 
contrast to simple parametric forms for dar k energy equation 
of sta te (EOS), such as w(z) - wp + w' z ("Cooray & Huterei| 
19991) and wC z) = wq+ WaZl{\ + z) dCheval lier & Polarski 200ll 



Lind er 2003) , a nearly model-independent approach was intro- 
duced by Huterer & Cooray (2005) to estimate the evolution of 
dark energy independently, and was adopted in analyses using 
type la supernova (SN la) data ("Riess et al."2007'; 'Sulliv an et alj 
[2007; Sarkar et al. 2008). Oi et al. (2008) extended the ap proac h 
to gamma-ray burst (GRB) luminosity data (iSchaefeiJ l2007h . 
for which higher redshifts are accessible compared to SNe la. 
Though more stringent results have been obtained by this, the 
dark energy EOS at high redshifts, where we only have GRBs, 
is still totally unconstrained, except that it is most likely nega- 
tive. This is primarily due to matter dominating dark energy at 
high redshifts and to dark energy becoming less important for 
determining the cosmic expansion; therefore, constraining it be- 
comes more difficult. In this work we explore both how many 
GRBs are needed to get substantial constraints on dark energy 
EOS at high redshifts beyond those of SNe la and where the 
GRBs' contributions lie for a foreseeable number of GRBs. 



2. Methodology 

For the uncorrected estimates of dark energy EOS, we first sep- 
arate the redshifts into several bins and assume a constant dark 
energy EOS for each bin. Then we generate a Markov chain 
through the Monte-Carlo algorithm according to the likelihood 
function. The covariance matrix of the dark energy EOS param- 
eters is calculated based on the Markov chain, and a transforma- 
tion derived from the covariance matrix to decorrelate the EOS 
parameters. The evolution of dark energy i s finally estimated 
out of the uncoiTelated EOS parameters (see lHuterer & Cooravl 
(120051) for details). 

The process of generating Markov chains in these proce- 
dures is very time-consuming, especially when there are con- 
siderable observational/mock standard candles. We see that we 
need a large number of GRBs to get substantial constraints on 
dark energy EOS at high redshifts. Instead of using the uncorre- 
cted approach directly, we adopt an alternative way to estimate 
this number. When estimating the constraints imposed on dark 
energy EOS in a certain redshift range, we allow the EOS param- 
eter to vary only in that redshift bin and fix the EOS parameters 
elsewhere, say, to -1. This is compatible with the uncoiTelated 
approach to an approximat e estimate, because the un coiTelated 
EOS parameters defined in iHuterer & Cooravl (l2005l) are them- 
selves localized well in redshift bins. Of course, the fixation of 
EOS at some redshifts will cause underestimates of the errors 
on the EOS at other redshifts, so the number we have estimated 
here can be viewed as a lower limit. 
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We determine the likelihood for EOS parameters from the;^f^ 
statistic. 



An 



!=1 



cr? 



(1) 



where jUo is the observational/mock distance modulus of the stan- 
dard candles with standard deviation cr, ^,{0) is the theoretically 
predicted distance modulus based on a cosmological model with 
parameter set 9, and n is the number of the standard candles. We 
constrain dark energy EOS parameters with mock data, whose 
generation involves random numbers (see the details below). To 
reduce the impact from the fluctuation of the mock data them- 
selves, we generate many more standard candles than needed. 
For example, if we want to see the constraints imposed by n sam- 
ples of standard candles, we actually generate A^ mock samples 
with N » n and calculate ;t'^. Then;^'^ is given by 



/vn M^^' 



(2) 



Throughout this paper, we choose A^ = 1.0E5. 

The fiducial cosmological model we used to generate mock 
data is the flat ACDM model with Q.„, = 0.279 and //q = 
70.1km/s/Mpc (see K omatsu et al. (2008) ). For SNe la, we ap- 
proximate the intrinsic noise in distance moduli as a Gaussian 
scatter with a dispersion of 0.1 magnitudes and the observa- 
tional errors as 0.23 magnitudes, which are approximately the 
average s of the observational errors of the presen t ly available 
SNe la (jRiess et alj|2007l: IWood-Vasev et al.1 120071; iDavis et alj 
|2007|) . The total error in our generated distance moduli for SN 
la is therefore VO.l^ + 0.23^ magnitudes. We assume an uni- 
form distribution for SNe la along the redshifts. For GRBs, in- 
stea d of generating mock data about the five luminosity relations 
(see lSchaefed (1200 7)). we directly generate distance moduli like 
we do for SNe la for simplicity. The intrinsic scatter is set to be 
0.65 magnitudes, which is approximately the average of the er- 
rors of the GRBs' average distance moduli presented in lSchaeferl 
(|2007|) . and we ignore the measurement uncertainties, which are 
less than the intrinsic scatter. We consider two kinds of distribu- 
tions for GRBs in the redshift bin 1.8 < z < 7. One is the uni- 
form distribution, the other a very rough approximation t o the 
distribution presented by Fig. 2 in iBromm & Loebl (l2006h . i.e. 
P{z) oc exp(-z/7). We will see that our results are independent 
of the GRB distributions. 
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Fig. 1. Constraints from GRBs uniformly distributed in the red- 
shift bin 1.8 < z < 7 on the dark energy EOS parameter 
w(1.8 < z < 1). The different lines stand for different number 
of GRBs. The probabilities are normalized to be 1 at the max- 
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Fig. 2. Constraints from GRBs distributed in the redshift bin 
1.8 < z < 7 according to P{z) oc exp(-z/7) on the dark energy 
EOS parameter w(1.8 <z<l). 



3. Results 

Figures [U and |2] show our results for the constraints from GRBs 
distributed in the redshift bin 1.8 < z < 7 on the dark energy 
EOS parameter w(1.8 < z < 7). We can see that, for a few hun- 
dred GRBs, the constraints are only the steep drop at about zero 
that is seen in the probability function of the EO S parameter 
P{w). This is consistent with the results in lOi et al] ( 12008.) . Only 
when we have more than about 5000 GRBs can we begin to get 
concrete constraints on the EOS parameter. The GRBs' distribu- 
tions have little impact on the conclusion; i.e., it is difficult to 
constrain the dark energy EOS parameters beyond the redshifts 
of SNe la with GRBs unless some new luminosity relations for 
GRBs with smaller scatters are discovered. 

However, this does not mean that high-redshift GRBs con- 
tribute little to constrainin g the dark energ y EOS parameters. It 
has been demonstrated in lOi et al.l ( 20081) tha t, even with the 
presently available 69 GRBs (ISchaeferl |2007|) . the constraints 



could be improved significantly at redshifts 0.5 < z < 1.8. 
Part of the improvement stems from GRBs beyond redshift 1.8. 
Because the luminosity distances of standard candles depend on 
the behavior of the dark energy through an integration over the 
redshift, high-redshift GRBs put constraints on dark energy at 
lower redshifts, where dark energy is important in determining 
the cosmic expansion. And since there are few GRBs at low red- 
shifts, the contributions from GRBs would lie primarily in the 
middle redshifts. In Figs. |3]and|4]we explicitly show the con- 
straints from GRBs distributed in the redshift bin 1 .8 < z < 7 
on the dark energy EOS parameter w(0.5 < z < 1.8). For com- 
parison, we also plot the constraints from SNe la uniformly dis- 
tributed in the redshift bin 0.5 < z < 1.8 on the dark energy 
EOS parameter vv'(0.5 < z < 1.8) in Fig.|5] We can see that the 
contributions from GRBs are comparable to that from SNe la. 
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4. Summary 
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Fig. 3. Constraints from GRBs uniformly distributed in the red- 
shift bin 1.8 < z < 7 on the dark energy EOS parameter 
w(0.5<z< 1.8). 
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Fig. 4. Constraints from GRBs distributed in the redshift bin 
1 .8 < z < 7 according to P{z) oc exp(-z/7) on the dark energy 
EOS parameter w(0.5 <z< 1.8). 
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We explored the GRBs' contributions in constraining the dark 
energy EOS at high redshifts (1.8 < z < 7) and at middle red- 
shifts (0.5 < z < 1.8). When constraining the dark energy EOS 
in a certain redshift range, we allow the dark energy EOS pa- 
rameter to vary only in that redshift bin and fix EOS param- 
eters elsewhere to -1. We find that it is difficult to constrain 
the dark energy EOS parameters beyond the redshifts of SNe 
la with GRBs unless some new luminosity relations for GRBs 
with smaller scatters are discovered. However, at middle red- 
shifts, GRBs have contributions comparable with SNe la in con- 
straining the dark energy EOS. 
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Fig. 5. Constraints from SNe la uniformly distributed in the 
redshift bin 0.5 < z < 1.8 on the dark energy EOS parameter 
w(0.5 <z< 1.8). 



